Optical Control of a Quantum Rotor 
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The possibility to coherently control a quantum rotor is investigated theoretically. The rotor is 
realized by an antiferromagnetic spin-1 Bose-Einstein condensate, trapped in the optical field of a 
Fabry-Perot resonator. By tuning the pumping field of the resonator, coherent control over the rotor 
is achieved. The technique is illustrated by the numerical simulation of a protocol that transforms 
the rotor's ground state into a squeezed state. The detection of the squeezed state via measurement 
of intensity-correlations of the cavity field is proposed. 
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Among all quantum mechanical degrees of freedom, 
the spin stands out as a fundamental property with no 
classical analogue. While many effects of spin in crystals 
are well understood and readily seen in the macroscopic 
world through the magnetism of solids, the theoretical 
study of spin models continues to be a topic of active re- 
search. The realization of spinor Bose-Einstein conden- 
sates (BECs) pQ , with their remarkable control of dimen- 
sionality, external potentials and two-body interactions, 
has significantly broadened the experimental toolkit for 
the study of such models and already lead to a wide range 
of rich new physics, see [2] for a recent review. 

Besides fundamental aspects, spinor BECs also hold 
promises for ultra-sensitive magnetometry with large spa- 
tial resolution [SJ. In this context, much attention has 
been given to spin-squeezing [4H6] , a consequence of non- 
linear spin-mixing dynamics. Some groups have also re- 
ported the observation of other quantum signatures in 
the population dynamics of such gases [342] ■ However, 
in order to use spinor BECs as a nonclassical resource, co- 
herent control over the quantum state has to be achieved, 
which is difficult due to the intrinsically nonlinear nature 
of the system. 

This letter proposes a scheme that allows coherent 
control over a quantum rotor, realized by an anti- 
ferromagnetic spin-1 BEC trapped within the field of an 
optical resonator [TQl - fT2] . We show how the optical field 
provides a convenient knob to coherently control the dy- 
namics. As an example, we demonstrate the possibility 
to squeeze the quantum rotor by temporally modulating 
the pump-field of the optical cavity. Our results arc more 
far-reaching however, as they also provide a new mech- 
anism to produce spin squeezing, since control over the 
rotor amounts to control over the spin degrees of freedom 
of the spinor gas. 

We consider the spin-mixing dynamics of an anti- 
ferromagnetic spin-1 BEC in the case of a tightly trapped 
cloud, where all three spin populations can be described 
by a single spatial mode. The spin degree of freedom ex- 
hibits several peculiarities that make it particularly suit- 
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FIG. 1: A spin-1 BEC is trapped in an optical cavity and 
dispersively coupled to the linearly polarized cavity mode. 
The mf = transition is forbidden by selection rules. 



able for the preparation of non-classical states. Firstly, 
the internal time scales of spin systems are very slow 
compared to optical time scales. Secondly, the spin of a 
BEC is very well isolated from decoherence and thermal 
effects, with decoherence times in the range of seconds. 
This isolation is a natural requirement to perform coher- 
ent manipulations on the system. 

Before including the light field in our treatment, let 
us focus on the condensate alone. Its spin populations 
evolve due to spin-dependent collisions and the linear 
and quadratic Zeeman shifts. Since the two shifts com- 
mute and evolve on vastly different timcscalcs, they can 
be studied independently. In our case, we focus on the 
quadratic shift alone. The Hamiltonian describing the 
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system is 

#coi - °^F 2 - gSJSo, (1) 

where c 2 is the spin coupling, N the total particle number 
and F = b^F^bj, with i 7 ^- spin-1 matrices and bj bosonic 
annihilation operators for the three m,f = j £ {±,0} 
components. The second term is the quadratic Zeeman 
shift due to an external magnetic field, with q being pro- 
portional to its squared absolute value. 

It was shown in Ref. [12] that the collision part of 
Hamiltonian (nl can be mapped onto an angular mo- 
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mentum Hamiltonian, c 2 1L /(2N), acting on the Eu- 
clidean space spanned by \x) = (|+) — |— })/V2, \y) = 
(1+) + h))/(*V2) and \z) = |0), where {|±),|0)> is 
the spin basis. Physically, the orientation of the ro- 
tor gives the axis in real coordinate space along which 
the z-component of the spin vanishes. More recently it 
was shown how to extend that mapping to include the 
quadratic Zeeman shift |101 111] , resulting in an external 
potential V{6) 

Hrot = \j+qV{6), (2a) 
V(§) = X ism 2 09)+ X2 sin 2 (2f?). (2b) 

Here X i = (N + 3/2), X 2 = qN/(8c 2 ) and I = Nh 2 /c 2 
and functions of operators are defined by their power 
series. The ratio between the two terms in the potential 
is given by X1/X2 ~ <?/(8 c 2) and measures the balance 
between the Zeeman shift and the collision term. For an 
external field around 100 mG and a sodium condensate, 
X1/X2 ~ 10 -3 , and we can neglect the term proportional 
to X2, rendering the V{&) proportional to sin 2 (#). 

As a result of the symmetry of (JlJ under exchange 
of mf — ±1 atoms, the angular momentum along the 
z-axis, which in terms of atomic variables is given by 
is a conserved quantity. Preparing 
the system initially with the same number of atoms in 
these two states reduces the number of degrees of freedom 
to only two. The two conjugate observables are given by 
the angle 8, measuring the population in the to/ = 
state, and the associated angular momentum Lg, related 
to coherences between the m j = and the other two spin 
components. Note that 9 is a periodic observable, but it 
is not related to the quantum phase |10j. 

We now show how it is possible to dynamically con- 
trol the state of the spinor gas and generate nonclassical 
rotor states by trapping it in an optical resonator driven 
by a laser of time-varying amplitude, see Fig. [I] If the 
laser field is linearly polarized along the direction of an 
external magnetic field and far detuned from the atomic 
transition, the Hamiltonian for the compound system in 



the frame rotating at the laser frequency oji reads [T2"] 

L 2 

fltot = HAa^a+-^+(q + MJ a) a) V(e)-ihT 1 (a-a^)+H K , 

(3) 

where the operator a annihilates a photon in the cav- 
ity mode and the atom-photon coupling is given by 
Uq = g 2 /(loi — uj a ) with g the dipole coupling and ui a the 
frequency of the atomic transition. We assume Uo > 
without loss of generality and also have introduced the 
pumping rate 77, the cavity detuning A = lo c — loi with 
lu c the resonance frequency of the empty cavity and H K 
accounting for optical dissipation with rate n. Dissipa- 
tion for the atomic system occurs on the timescale of 
seconds, while the dynamics considered here take place 
within milliseconds and we can safely neglect atomic dis- 
sipation. The evolution of the angular operator 6 is given 
by 

d 2 9 

I-jjp = -2(q + hU a a^a) 

x (xi sin(0) cos(§) + 2 X 2 sin(2(9) cos(2(9)) . (4) 

In the following we treat the driving field classically, and 
due to the slow evolution of the spin dynamics we can 
assume the optical field to follow the evolution of 6 adia- 
batically. Eliminating the optical field dynamics results 
in the effective potential 

M 0- + ( j*±jwg& Y (5) 

qK \ K I 

and the effective Hamiltonian for the one-dimensional ro- 
tor becomes 

n eS = ^+qV eS (9). (6) 

Neglecting quantum fluctuations in the cavity field is 
not a priori unproblematic since for practical parameters 
the average number of photons in the cavity can be rather 
small, e.g. less than a thousand. The effective potential 
still provides a correct description for the rotor because 
the fluctuations in the cavity are correlated on a timescale 
of k^ 1 , much faster than the characteristic timescales of 
the rotor. The rotor averages out quantum fluctuations. 
This is confirmed in our numerical calculations, where we 
have included fluctuations of the optical field. Before we 
present these results, however, let us discuss the merits of 
the optical field as a control tool to coherently manipulate 
the rotor. 

The qualitative nature of V c s depends on the rela- 
tions between A, UqN and k. The inverse tangent is 
a monotonously increasing function that changes very 
slowly, except in an interval of order unity around the 
origin, where it increases steeply. Thus the minima and 
maxima of the original potential remain conserved. We 
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FIG. 2: Effective scaled rotor potential as a function of cavity- 
detuning and angle. See text for list of parameters. 

can moreover expect that if the absolute value of the 
argument is very large, the qualitative distortion of the 
original potential is small. In terms of our parameters, 
such a regime is captured by the condition |A| ^> UqN , 
corresponding to the situation of a pump laser far de- 
tuned from the cavity resonance for any state of the ro- 
tor. In this case we can expand the potential around 
Uq(V(0)) I k and find after dropping constant terms 

which is the trapping enhancement found in [12] . 

For the resonator to significantly distort rather than 
merely scale the potential, two conditions need to be met: 

1. The total range of detuning due to the rotors state, 
UqN / n, has to be much larger than unity. 

2. The argument of the inverse tangent has to cover 
regions around the origin, where the function 
changes characteristics. 

The first condition requires a sufficiently strong atom- 
photon coupling, while the second condition can be sat- 
isfied by the right choice of laser frequency, i.e. A < n 
and |A| < UqN. If both conditions are satisfied, the op- 
tical resonance will enhance the parts of V(9) that fall 
on the steep portion of the inverse tangent around zero 
and flatten those beyond this region. For A = 0, the 
rotor's confinement is maximally tightened, whereas for 
A = — UoN we are left with a predominantly flat ef- 
fective potential. This behavior is illustrated in Fig. [2] 
which shows the scaled potential as a function of 9 and 
the pump-cavity detuning A. 

The character of the effective potential will adapt to 
changes in 77 and A on the fast time scale K , but we are 
free to change them as slowly as we want. This flexibil- 
ity allows for efficient manipulations of the rotor's quan- 
tum state. For the remainder of this discussion, we will 



turn to the Schrodinger picture in space, which we also 
used for our numerical calculations. Let us assume for 
concreteness that the quantum rotor is initially in some 
state 

iv>w) = E c ^'W>' (8) 
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where \4>j{&)) are the eigenfunctions corresponding to a 
particular form of the effective potential. Under a change 
of V e ff sufficiently slow that the adiabatic theorem holds, 
the coefficients Cj do not change their distribution, but 
the eigenfunctions tpj{9) become eigenfunctions ipj(9) °f 
the new potential. This can be used to reset the system to 
its original parameters or to bring it to a desired regime of 
operation. For a diabatic change, however, the new state 
is the projection of the old one on the new eigenfunctions, 

wm=Y f c j wmHo)Wi)- w 

Immediately following the switch the state still has the 
same wave function, but with a different distribution of 
eigenfunctions. Diabatic changes redistribute the weight 
of the eigenvalues of a wave function and hence the state's 
characteristics. 

With appropriate combinations of adiabatic and dia- 
batic transitions it is possible to engineer desirable states, 
as we now illustrate on the concrete example of squeezed 
state generation. The idea is to diabatically switch be- 
tween tight and wide potentials at the right instants in 
time and let the rotor evolve freely between the changes. 
This protocol was originally proposed to squeeze the mo- 
tion of trapped ions |14j . but the fast switching times 
required make it challenging in such systems |15j . For 
typical experimental parameters the angular width of the 
ground state even in the shallow potential is very small, 
so that we can neglect the periodicity of the rotor and ap- 
proximate both regimes by harmonic confinements with 
frequencies uj\ and oj 2 for the tight and wide trap respec- 
tively - these being easily found by an expansion of the 
potential. 

We assume that the rotor is initially in its ground state 
with respect to the trapping frequency oji. The distri- 
butions of 9 and Lg are Gaussian with their respective 
zero-point widths. Suddenly changing the rotor's trap- 
ping frequency to the lower frequency uj^ leaves the wave- 
function unchanged; with respect to the new potential, 
however, its angular distribution is narrower than the 
zero-point width while it's angular momentum distribu- 
tion is wider, i.e. the rotor is squeezed with squeezing 
angle in the — Lg plane. In the new potential, the 
rotor undergoes its free time evolution and after a time 
t = 7r/2o;2 the angle is converted to angular momentum 
and vice versa. Now switching back diabatically to the 
initial, tight potential the situation is exactly reversed 
compared to the first change and the rotor experiences 



further squeezing, the squeezing angle now being tt/2. 
After the time t — n/2uii, this angle has changed to ir 
and we again switch to the wide potential. This process 
can be repeated and after each cycle, the uncertainty in 
one quadrature is scaled by the factor (cj 2 /wi) 2 , while the 
orthogonal quadrature is scaled with the inverse factor. 
This protocol can be understood as a discretized version 
of the squeezing Hamiltonian a 2 + (a^) 2 that is realized 
for example in parametric amplifiers. 

In practice, things differ of course from this idealized 
picture. First, it is impossible to perform truly instanta- 
neous changes to the potential. Moreover, the timing of 
the switches is essential, but since typically the number 
of particles in the condensed gas is not constant dur- 
ing an experiment and also exhibits shot-to-shot fluctu- 
ations, the rotation frequencies of the rotor are not pre- 
cisely known. A further limiting point is the anharmonic 
nature of the trapping potentials, which have a sin 2 (#) 
dependence. These last two points become more relevant 
as more squeezing is achieved and the rotor's distribution 
becomes wider. 

To investigate these difficulties we have performed a 
numerical simulation of the protocol. Specifically, we 
solved the time dependent Schrodinger equation in 6- 
space with a realistic time dependence of the effective 
potential, including the effects of photon noise through 
stochastic fluctuations of the effective potential and also 
chosen the atom number with a 5% uncertainty to mimic 
potential experimental difficulties. The result of these 
simulations is shown in Fig. [3]for a 23 Na BEC consisting 
of N « 10 4 atoms with a number density of 5 x 10 14 cm~ 3 . 
The pump laser is detuned from the sodium Dl line by 
A a = 2n x 1 GHz and feeds a cavity with line-width 
K = 2tt x 1 MHz. Using these parameters, we find that 
the maximal confinement and the flattest effective poten- 
tials have harmonic frequencies of u>i ~ 2ir x 43kHz and 
cj2 ~ 2n x 7kHz respectively, separating the rotor and 
cavity timescales by more than a factor of 20. The rotor 
was initially in its ground state, see panel (a), and the 
squeezing after the potential switches is clearly visible. 

There are several possibilities for the experimental 
characterization of the state of the rotor. The most 
straightforward approach is to perform a Stern-Gerlach 
type measurement, which is a projection on the number 
basis of the three components. While the angular dis- 
tribution is directly related to the population of m/ =0 
atoms, the rotor's angular momentum corresponds to co- 
herences between mf = and rrif = ±1 atoms and thus 
is not captures in this measurement. Furthermore, it is 
a destructive measurement technique, and thus not suit- 
able to acquire large amounts of data. To fully resolve 
the rotor dynamics could involve dispersive imaging tech- 
niques that would give direct in-situ access to atomic co- 
herences [TB]. Continuous monitoring of the coherences, 
however, will generally require quantifying the measure- 
ment back-action on the rotor |17j . Moreover, additional 
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FIG. 3: Wigner distributions of the rotor, (a) Ground-state 
distribution, (b) Distribution after the first switch of the 
trapping potential, (c) Distribution after completion of the 
squeezing cycle. See text body for a list of system parameters. 
Units are the zero-point widths of the system. 



imaging beams impose experimental difficulties. 

The necessary information is also contained in the out- 
put field of the optical resonator, provided that one ac- 
counts for its quantum nature. Assuming that the ro- 
tor does not change its state over the fast time scale 
we can expand the potential in Hamiltonian ([3| as 
V(0) = (V(0)) + SV and solve the equation of motion for 
the photon field operator to first order in 5V 



a(t) =a s + F in + jSV, 



(10) 



where a s is the coherent steady-state amplitude, F- m in- 
volves the vacuum input noise and 



7 = ir lUo 



K/2 + i{A + U Q (V(9))) 



(11) 



Explicit expressions are easily found from Eq. ^ , but are 
omitted here for brevity. From these expressions one can 
find a non- vanishing contribution of (<5U 2 ) to the second 
order correlation function </ 2 ) = (aJaJaa)/[{aJa){aJa)]. 
Our numerical calculations indicate that after 5 squeez- 
ing cycles, which would take about 200 [isec, the state of 
the rotor would lead to oscillations of with an am- 
plitude of order ~ 0.02, which would render our proposal 
realizable in state-of-the-art experiments. 

In summary, we have investigated the possibility to 
control the quantum state of a spinor gas through an 
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optical resonator. Using the rotor mapping for the quan- 
tum gas, we find the effective potential created by the 
optical field. Due to the slow spin-mixing dynamics, dia- 
batic changes in the effective potential can be induced by 
tuning optical parameters. We demonstrated how such 
changes can be used to squeeze the rotor by means of 
a numerical simulation. Information about the rotor's 
state is imprinted on the intensity fluctuations of the 
cavity field and may be used to verify squeezing. Further 
studies will involve the influence of in-situ measurements 
on the rotor, as well as the possibility to couple separate 
spinor clouds through the cavity field. 
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